Running of a s in the MSSM with three-loop accuracy 
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Abstract. The evolution of the strong coupling constant a s from Mz to the GUT scale is presented, 
involving three-loop running and two-loop decoupling. Accordingly, the two-loop transition from 
the MS to the DR scheme is properly taken into account. We find that the three- loop effects are 
comparable to the experimental uncertainty for a s . 

PACS. 11.30.Pb Supersymmetry - 12.38.-t Quantum chromodynamics 



1 Introduction 

Supersymmetry (SUSY) is currently believed to play 
an important role in physics beyond the Standard Model 
(SM). There are several reasons that point out the 
Minimal Supersymmetric Standard Model(MSSM) as 
a preferred theory describing new physics. The fore- 
most ones are its milder divergency structure that solves 
the naturalness problem and the possibility to explain 
the electroweak supersymmetry breaking as a conse- 
quence of radiative corrections. Another compelling 
argument in favour of SUSY is the particle content 
of the MSSM, that leads in a natural way to the uni- 
fication of the three gauge couplings at a high energy 
scale fj, ~ 10 16 GeV, in agreement with Grand Unifica- 
tion Theories (GUT). This observation together with 
the consistent predictions made for SM parameters, 
such as the top quark mass and the ratio of the bot- 
tom quark to the tau lepton masses, using constraints 
on the Yukawa sector of SUSY-GUT models, brought 
SUSY in the center of the phenomenological studies. 

Nevertheless, SUSY can only be an approximate 
symmetry in nature and several scenarios for the mech- 
anism of SUSY breaking have been proposed. A possi- 
bility to constrain the type and scale of SUSY break- 
ing is to study, with very high precision, the relations 
between the MSSM parameters evaluated at the elec- 
troweak and the GUT scales. The extrapolation over 
many orders of magnitude requires high-precision ex- 
perimental data at the low energy scale. A first set of 
precision measurements is expected from the CERN 
Large Hadron Collider (LHC) with an accuracy at the 
percent level. A comprehensive high-precision analysis 
can be performed at the International Linear Collider 
(ILC), for which the estimated experimental accuracy 
is at the per mill level. In this respect, it is neces- 
sary that the same precision is reached also on the 
theory side in order to match with the data p] . Run- 
ning analyses based on full two-loop renormalization 
group equations (RGEs) [2H3] for all parameters and 



one-loop threshold corrections [1] are currently imple- 
mented in the public programs ISA JET [5], SOFT- 
SUSY [6j, SPHENO [7], SuSpect 0. The agreement 
between the different codes is in general within one 
percent [9j. A first three- loop running analysis, based, 
however, only on one-loop threshold effects, was car- 
ried out in Ref. [TP] . 

In this talk, we report on the evaluation of the 
strong coupling a s in MSSM, based on three-loop RGEs 
[TT] and two-loop threshold corrections [12j . On the 
one hand, the three-loop corrections reduce signifi- 
cantly the dependence on the scale at which heavy 
particles are integrated out [13]. On the other hand, 
they are essential for phenomenological studies, be- 
cause they are as large as, or greater than, the ef- 
fects induced by the current experimental accuracy of 
a s (Mz) [H]. Additionally, we compare the predictions 
obtained within the above mentioned approach with 
those based on the leading-logarithmic (LL) approxi- 
mation suggested in Ref. [T]. 



2 Evaluation of q; s (//gut) from a s (M z ) 

The aim of this study is to compute a s at a high- 
energy scale /i ~ O(^gut), starting from the strong 
coupling constant at the mass of the Z boson Mz- 

We denote this parameter af S '^ 5 \Mz) to specify that 
the underlying theory is QCD with five active flavours 
and MS is the renormalization scheme. The value of 
&s(hgut) is the strong coupling constant in the MSSM 
renormalized in the DR-scheme, that we denote as 

DR.(full)/ n „ , r DR. (full) , 

a s (a*gut)- lor the evaluation ot a s from 

a MS,(n f ) we f n ow the "common scale approach" [T5] . 
which requires a unique scale for the matching between 
QCD and MSSM. More precisely, for mass indepen- 
dent renormalization schemes like MS or DR, the de- 
coupling of heavy particles has to be performed ex- 
plicitely. In practice, this means that intermediate ef- 
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fective theories are introduced by integrating out the 
heavy degrees of freedom. One may separately inte- 
grate out every particle at its individual threshold ( "step 
approximation"), a method suited for SUSY models 
with a severely split mass spectrum. But the interme- 
diate effective theories with "smaller" symmetry raise 
the problem of introducing new couplings, each gov- 
erned by its own RGE. To overcome this difficulty, for 
SUSY models with roughly degenerate mass spectrum 
at the scale M, one can consider the MSSM as the full 
theory that is valid from the GUT scale fiGUT down 
to M, which we assume to be around f TeV. 

For the running analysis of the strong coupling con- 
stant, we can distinguish four individual steps that we 
detail below. 

1. Running of a^ IS '^" / ' ) from /j = Mz to fi = /idec- 
The energy dependence of the strong coupling con- 
stant is governed by the RGEs 



be done at an arbitrary decoupling scale \x: 



/?(a s ) = -a^/3„cC 

n>0 



(1) 



In QCD with n/ quark flavours, the f3 function is 
known through four loops both in the MS [T6UT7] 
and the DR-scheme |18j . 

„ „ r MS,(j»f ) ; n , DR.Of), n 

2. Conversion of a s (A*dec) to a s (/^decj- 
For the three-loop running analysis we are focusing 
on, one needs to evaluate the dependence of a s 
values in the DR scheme from those evaluated in 
MS scheme with two loops accuracy [18] 



MS 



,DR 



V DR 



4tt 
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(2) 



Here, the following notations have been used a 



DR. 



a? R ' (n/) (M) anda^=- MS ' ( " /) 



(fi). a e = ae (//) 



is one of the so-called evanescent coupling con- 
stants that occur when DR is applied to non super- 
symmetric theories (QCD in this case). In particu- 
lar, it describes the coupling of the 2e-dimensional 
components (so-called £-scalars) of the gluon to 
a quark. It is an unphysical parameter that must 
decouple from any prediction for physical observ- 
ables. We also used this property as a consistency 
check for our method. As mentioned above, we as- 
sume that QCD is obtained by integrating out the 
heavy degrees of freedom (squarks and gluinos) 
from SUSY-QCD. In this case, the evanescent cou- 
plings are uniquely determined by the matching 
conditions between the two theories, that we dis- 
cuss below. 



3. Matching of a? R '' n/ ^(^dc 



, DR,(full), n 

and a s (/idee ) 



Integrating out all SUSY particles at the common 
scale of SUSY mass spectrum, one directly obtains 
the SM as the effective theory, valid at low ener- 
gies. The transition between the two theories can 



a 



DR.,(™/) 



(/*) = 



a 



DR,(full) 



(/<) 



(3) 



and £ e depend logarithmically on the scale fx, 
which is why one generally chooses /x ~ M. In 
Eq. ©, nj = 6 means that only the SUSY par- 
ticles are integrated out, while for nj — 5 at the 
same time the top quark is integrated out. 
In a supersymmetric theory, SUSY requires that 

a.e nll \n) — a? R '^ ful ^(/i) at any scale. Let us re- 
mark that a^ R " (//dec) is not known a priori 
and one cannot use Eq. ([3]) directly in order to 

derive a e ■ Rather, we start with a trial value 
for a? R '^" H ' (/idec) and obtain the corresponding 
a e 5 ^(/Xdcc) as weu as Qi? (Mdec) through Eq. (|3|). 

MS, 

MS ' {5) (M Z ). The trial value for 



Then we evaluate Ohf (fidec) through Eq 
and from that a 

DR, (full) 



(Udec) is systematically varied until the 

resulting a^'^^M^) agrees with the experimen- 
tal input. 

4. Running of a]j )R " fu1 ^ from /i = /idee to /i = /iQUT- 
The energy dependence of the strong coupling con- 
stant is in this case governed by the MSSM RGEs. 
In SUSY-QCD, the f3 function has been evaluated 
in the DR-scheme through three loops [IT] . 

Assembling the above mentioned steps, we can pre- 
dict the value of a s (/icuT) with up to three- loop ac- 
curacy. This procedure is implemented in most of the 
present codes computing the SUSY spectrum [TUSHE] 
by applying the one-loop approximation of Eq. ([3|) and 
setting rif = 5 and fi = Mz- The advantage of this pro- 
cedure as compared to a multi-scale approach is that 
the RGEs are only one-dimensional and that for a e 
one can apply Eq. ©. 

Let us note that in principle it is possible to decou- 
ple the top quark separately. The only new ingredient 
needed is the decoupling constant for going from five 
to six quark flavours in the MS scheme. In any case, 
for a mass spectrum as given by the benchmark point 
SPSla' pQ, for example, the separate decoupling of the 
top quark implies a numerically small effect. This can 
also be established by comparing "Scenario D" and 
"Scenario C" in Ref. [12]. 

The phenomenological significance of the three-loop 
order corrections is discussed in detail in the next sec- 
tion. 



3 Numerical results 

The result for a° R (/z G UT = 10 16 GeV), obtained using 
M z = 91.1876 GeV, m t = 170.9±1.9 GeV , ap(M z ) = 
0.1189 , and M = m$ = m§ = 1000 GeV as input pa- 
rameters is shown if Figure [1] The dotted, dashed and 
solid line are based on the approach described above, 
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Fig. 2. as(uGUT) as a function of M. 



non-unification behaviour of the three couplings at the 
GUT scale. 



4 Conclusions 

We have used recent three- and two-loop results for 
the (3 functions and the decoupling coefficients, respec- 
tively, in order to derive a° R (/iGUT) from af^M^) 
at three-loop level. 

It turns out that the three-loop terms are numeri- 
cally significant. The dependence on where the SUSY 
spectrum is decoupled becomes particularly flat in this 
case. The theoretical uncertainty is expected to be neg- 
ligible w.r.t. the uncertainty induced by the experi- 
mental input values. 

Comparing our results and methods to the liter- 
ature, we find that the issue of evanescent couplings 
has either been ignored (by assuming a e = a s ) or 
circumvented by decoupling the SUSY spectrum at 
Hdcc = Mz- We find that at one- and two- loop level, 
this choice does not allow for a good approximation 
of the higher order effects, if one assumes the SUSY 
partner masses to be of the order of 1 TeV. 

In consequence, we recommend that phenomeno- 
logical studies concerning the implications of low en- 
ergy data on Grand Unification should be done at 
three-loop level. 
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where n-loop running is combined with [n— l)-loop de- 
coupling, as it is required for consistency (n = 1,2,3, 
respectively). We find a nice convergence when going 
from one to three loops, with a very weakly /idee - 
dependent result at three-loop order. For comparison, 
we show the result (the dash-dotted line) obtained 
from the formula given in Eq. (21) of Ref pQ. It corre- 
sponds to the resummed one-loop contributions origi- 
nating from both the change of scheme and the decou- 
pling of heavy particles. However, the difference be- 
tween our three-loop result with two-loop decoupling 
(upper solid line) and the one-loop formula given in 
Ref. p] exceeds the experimental uncertainty by al- 
most a factor of four for sensible values of /idee- This 
uncertainty is indicated by the hatched band, derived 
from 8a s {M z ) = ±0.001 [HJ. The formulae of Ref. [TJ 
should therefore be taken only as rough estimates. 

In Figure [2] we show Q! s (//gut) as a function of M 
where /^dec = M has been adopted. Dotted, dashed 
and full curve correspond again to the one-, two- and 
three-loop analysis and the uncertainty form a s {Mz) 
is indicated by the hatched band. One observes a vari- 
ation of 10% as M is varied between 100 GeV and 
10 TeV. This shows that the actual SUSY scale can 
significantly influence the unification, respectively, the 
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